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ABSTRACT 

Aims. We study the spectroscopic binary system Gl 375. 

Methods. We employ medium resolution echelle spectra obtained at the 2.15 m telescope at the Argentinian observatory 
CASLEO and photometric observations obtained from the ASAS database. 

Results. We separate the composite spectra into those corresponding to both components. The separated spectra allow 
us to confirm that the spectral types of both components are similar (dMe3.5) and to obtain precise measurements 
of the orbital period (P = 1.87844 days), minimum masses (Mi sin^ i = 0.35 M0 and M2 sin^ i = 0.33 Mq) and other 
orbital parameters. The photometric observations exhibit a sinusoidal variation with the same period as the orbital 
period. We interpret this as signs of active regions carried along with rotation in a tidally synchronised system, and 
study the evolution of the amplitude of the modulation in longer timescales. Together with the mean magnitude, the 
modulation exhibits a roughly cyclic variation with a period of around 800 days. This periodicity is also found in the 
flux of the Ca II K lines of both components, which seem to be in phase. 

Conclusions. The periodic changes in the three observables are interpreted as a sign of a stellar activity cycle. Both 
components appear to be in phase, which implies that they are magnetically connected. The measured cycle of ~ 2.2 
years (~ 800 days) is consistent with previous determinations of activity cycles in similar stars. 

Key words. <binaries: spectroscopic - Stars: activity - Stars: chromospheres - Stars: flare - Stars: fundamental parameters 
- Stars: late-type> 



1. Introduction 

Gl 37 5 (V* LU Vel, HIP 48 904) is a V = 11.27, B-V = 
1.56 (|Perrvman et al.1 |1997[ ) spectroscopic binary system 
located 16 pc away in the Southern constellation Vela, at 
RA (J2000) = 09:58:34 and DEC (J2000) = -46:25:30. Both 
components of the system exhibit large levels of chromo- 
spheric emission. 

Until recently, both the duplicity of the system and 
its high level of chromospheric activity were unknown. A 
considerable number of photometric measurements exist in 
the literature, including very early works (see for exam- 



tne uterature, inciuamg very early worKs (see tor exam- 
ple iKron et al.l.ll957t ICorben et atl . 119721: iDovle fc ButlMl . 



I1990I: iBessell |1990D . However, the flare activit y of Gl 375 



was d iscovered in a relatively recent work by lOovle et al.l 
(|l990t ) using optical spectroscopy and optical and infrared 
photometry. They reported a high rate of flare occurrence 
and classified Gl 375 as d M3.5e , similar to the spectral 
classification of 'Bid elmanI (jl985l ). who, however, did not 
indicate the presence of em ission lines. The spectral reso- 
lution of iDovle et al.l ()1990[ ) was not enough to resolve the 
double features in the spectra of Gl 375, but they noted 
that the FWHM of the Ha line was twice that found in 



* The authors are visiting astronomers of the Complejo 
Astronomico El Leoncito operated under agreement between the 
Consejo Nacional de Investigaciones Cientfficas y Tecnicas de la 
Repiiblica Argentina and the National Universities of La Plata, 
Cordoba and San Juan. 



the spectra of a similar star. Therefore, they concluded 
that G l 375 was either a fast rota t or or a binary. Only re- 
cently, IChristian fc Mathioudakid (|2002f) obtained a high- 
resolution optical and IR spectrum of Gl 375, which allowed 
them to identify the characteristic double features and re- 
port the binary nature of this object. 

Stars in close binary systems tend to show higher 
levels of activity than single stars with the sam e char- 
acteristics (see, for ex ample ISchriiver fc Zwaanl Il99ll : 
IZaaarashvili et all . 120021 ). Due to its vicinity and relative 
brightness, the present system constitutes an interesting 
opportunity to study two very active stars in interaction. 

The afl dynamo is accepted to be the cause of activ- 
ity cycles in late-type stars with an outer convective layer. 
This dynamo is thought to be the result of the action of 
differential rotation at the interface between the convective 
envelope and the radiative core. Therefore, the presence 
and characteristics of activity cycles are closely related to 
the existence and depth of an outer convection zone. Since 
this depth depends on spectral type -from F stars that 
have shallow convection zones to middle M stars that are 
totally convective-, it is of special interest to study these 
cycles in stars of different spectral types, and in particular 
in middle-M stars, to determine whether there is an onset 
of cyclic activity. Recently, an activity cycle has been de- 
tected in Proxima Centauri (jCincunegui et al.l . l2007a[ ). To 
our knowledge, this is the only fully convective M-type star 
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for which an activity cycle has been detected so far. This 
poses a question on the origin of such a cycle. 

The recent discovery of a 5 Earth- mass planet in the 
habitable zone of the M3 dwarf Gl 581 (jUdrv et all . I2007D 
has made the study of habitab ility in these stars increas- 
ingly relevant ( s ee, fo r example, iBuccino et al.1 . l2006l l2007t 
Ivon Bloh et al.l . [200l . Chan a;es in the rate of flare occur- 
rences and CMEs characteristic of activity cycles could in 
principle affect the habitability of a planet in close orbit. 

Gl 375 was observed as part of an ongoing program to 
study chromospheric activity in late-type main sequence 
stars. By the time of our first observations in the year 2000, 
its binary nature was unknown. However, it was quickly 
realized that this was a double system, with chromospheric 
emission which was similarly strong in both components. 

In this first paper we present medium-resolution spec- 
tra of Gl 375 spanning over six years. We obtain precise 
measurements of the orbital parameters of the system and 
find evidence for the presence of an activity cycle in both 
components. In Sect. 2 we describe the spectroscopic ob- 
servations and outline the reduction procedure. In Sect. 3 
we describe the method used to separate the spectra of 
each component and we present precise measurements of 
the orbital parameters of the system. In Sect. 4 we analyze 
photometric and spectroscopic data of Gl 375 and find pe- 
riodic variations consistent with an acticity cycle. Finally, 
in Sect. 5 we summarize our results and outline our conclu- 
sions. 

2. Observations and data reduction 

The spectroscopic observations were obtained using the 
REOSC spectrograph in the Jorge Sahade 2.15 m telescope 
at the Complejo Astronomico El Leoncito (CASLEO), lo- 
cated 2550 m above sea level in the Argentinian Andes. The 
spectra cover the entire optical range (from w 3800 to 6800 
A) in 24 echelle orders. We employed a 300 micron-width 
slit which provided a resolving power of i? = A/5A « 13200. 
The detector used was a 1024 x 1024 pixel TEK CCD cam- 
era. 

Observations began in 1999 and are still ongoing. We ob- 
tain about three observations a year, weather permitting. A 
log of the observations analyzed here is presented in Table[T] 
Each log entry actually consists of two individual observa- 
tions which are combined to remove cosmic rays, and the 
reported exposure time is the sum of the exposure times of 
the two individual observations. The individual spectra are 
used in Sect. 14.31 to look for flares in the observations. 

The observations were reduced and analyzed using stan- 
dard IRAfQ routines. We calibrate in flux the spectra, using 
low-dispersion long-slit spectra of the same star. For details 
of the re duction and flux calib r ation process we refer the 
reader to ICincunegui Sz MauasI (|2004D . 

The resulting spectral resolution is high enough to 
clearly resolve the double features in the spectra of Gl 375 
at appropriate phases. The exposure times ranged from 20 
minutes to 240 minutes, producing a wide range of S/N. 

Our spectral coverage allows us to study the effect 
of chromospheric activity in the whole optical spectrum 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



Table 1. Observation log for Gl 375 
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" In ddmmyy format 
' HJD-2450000 
" See Sect. 



simultaneously. In previous works we have studied the 
behaviou r of different line s in ou r complete stellar sam- 
ple (see ICincunegui et al.1 l2007bl for the H & K Ca II 
lin es and H^ , 'Di'az e t al.ri2007l for the Na I D lines, and 
.Buccino k Mauas . .2007f for the Mg 11 h and k UV lines). 

3. Spectra separation and orbital parameters 

The composite spectra were separated using the itera- 
tive method presented bv lGonzalez fc Levatol ()2006h . This 
method allows to compute the individual spectra and the 
radial velocities (RVs) of the two stellar components of the 
binary system. In each iteration the computed spectrum of 
one component is used to remove its spectral features from 
the observed spectra. The resulting single-lined spectra are 
then used to measure the RV of the remaining component 
and to compute its spectrum by combining them appropri- 
ately. 

We used as a template for the cross-correlations a spec- 
trum of a star of similar spectral type (AD Leo = Gl 388) 
observed with the same instrumental setup during the night 
of March 7 2004. The zero-point was established with an 
uncertainty of about 0.5 kms~^ by means of the observa- 
tion of RV standard stars. In all cases for RV measurements 
we excluded the emission lines. 

A mean spectra of very high S/N was obtained for each 
component, which were used to check the spectral type of 
the stars. These spectra are shown in Fig. [1] together with 
the spectrum of AD Leo we employed. To measure the spec- 
tral typ e of the components we employed the TiO bands 
used bv lPettersen fc HawlevI (jl989) to compute the intrin- 
sic brightness of a sample of red flare stars. Unfortunately, 
only three of these bands fall within our spectral range 
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Fig. 1. Mean spectra of the two components of the binary system Gl 375, together with the spectrum of AD Leo used 
as template. The entire spectral range observed is shown, and the spectra are displaced arbitrarily for clarity. As can be 
seen, the three spectra are remarkably similar, indicating that both the primary and secondary components of Gl 375 
have spectral type M3.5Ve. 



(4760 A, 4950 A, 5450 A). We carefully compared the 
strength of these bands in the mean spectra of each com- 
ponent and in AD Leo and found that both stars are very 
similar between the m, and similar to AD Leg , one of the 
dM3.5 stars used bv lPettersen fc HawlevI (|l989t ). Therefore, 
the spectral types of the components must be around M3.5. 
Additionally, we computed the absolute magnitude of each 
comp onent using the parall ax from the Hipparcos cata- 
logue (jPerrvman et al.l . [l997[ ) and our result that both stars 
are identical. We found that both components have abso- 
lute magnitude My = 11.01, which is in exce llent agreement 
with the value for AD L eo My = 10.98 fsee lJohnsoiill 19651: 
Ivan Altena et al.l . ll995t) . Unfortunately, we could not verify 
the obtai ned spectra l type direc tly by using the re l ations 
given by iReid et al.l (| 19951 ) and iKirkpatrick et al.l (|1991[) 
since the TiO bands employed in those papers lay outside 
our spectral range. 

Also evident from Fig. [T] is the high level of emission 
in the Balmer lines and in the H & K Ca II lines, which 
implies a high level of chromospheric activity. In Fig. [2] we 
compare the flux intensity of AD Leo and Gl 375 in Hq,, H^, 
and in the Ca II H & K lines. Note that the line intensities 
in the mean spectra of Gl 375A and Gl 375B are compa- 
rable to those from AD Leo, which is considered a very 
active flare star. This should be expected since, as men- 
tioned in the Introduction, stars in close binary systems 
usually exhibit enhanced levels of chromospheric activity 
when compared with single stars of the same character- 
istics (jSchrijver fc Zwaan, 1991; Zaqarashvili et al., 2002). 
In Sect. I4.3l we present evidence of a continuous interaction 
between the magnetospheres of the components which may 
lead to an enhanced level of chromospheric activity. 

Using the obtained values of RV we fit a Keplerian orbit 
to the system. In Fig. [3] we show the RV curves along with 
the resulting fit. The orbital parameters we obtained are 
presented in Tabled Note that since the minimum mass of 
the components is close to the expected value for M-type 
dwarves, which is around 0.4 Mq for spectral type M2 and 
0.22 Mo for M5 (|Coxl . l2000( l. the orbital inclination angle 



cannot differ much from 90°. Therefore the minimum dis- 
tance asini = (5.665 ± 0.035) R© should be similar to the 
actual separation between the components. For such a tight 
orbit, we wo uld expect to find the stars tidally locked. In 
this respect, IZahnI (jl977l ) studied the effects of turbulent 
viscosity in binary star systems and provided expressions 
for the characteristic time for orbital and rotational syn- 
chronization and orbital circularization. Employing Eqs. 6.1 
and 6.2 from that work, we find that for a system as the 
one studied here, the synchronisation time is 

W„c«(g)"' (a/i?i)6«8.9Myr 

and the circularization time is 

tcirc « (g(l + q)/2)-^ (a/i?i)' ~ 1.7 Gyr , 

where q — M2/M1 and Ri is the radius of the primary. 
For these calculat ions we hav e assumed the primary radius 
to be 0.4 i?0 fsee lCoxl[2000l ) and the orbital separation to 
be similar to its minimum value, a « 5.66 i?©. Therefore, 
since we found a value of the eccentricity which is undis- 
tinguishable from zero, we believe the system must have its 
orbital and rotational periods synchronised. However, to 
our knowledge no measurements of the period of rotation 
of this system exist in the literature. In Sect. l4.2l we provide 
further evidence that the rotational and orbital periods are 
indeed synchronised. 

4. Variability 

In this Section, we study the variability of Gl 375 as a 
whole and of each component. Besides using the spectra 
obtained at CASLEO, we also employed photometric obser- 
vations pro vided by the All Sky Automated Survey (ASAS, 
lPoimanskil . l2002f ). The ASAS data for this system cover the 
period between 2002 and 2006, with gaps during the year 
2002 and 2005. We discarded all the observations which 
were not qualified as either A or B in the ASAS database 
(i.e. we retained only the best quality data), as well as 7 out- 
lier observations. The final dataset consists of 295 points, 
with typical errors of around 15 mmag. 
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Fig. 2. Comparison of the Ha, H/3, Ca II H and Ca II K line emissions in AD Leo (solid line), Gl 375A (dashed line) and 
Gl 375B (dotted line). The spectra of Gl 375 are mean spectra obtained from the entire set of observations. 




Phase 

Fig. 3. RVs measurements phased to the orbital period P — 1.8784425 days. The squares represent measurements of the 
primary star and the circles are measurements of its companion. The solid line is the best fit to the data. In all cases the 
error bars are smaller than the symbols. In the lower panel we show the residuals of the fit and the individual errors. 



4.1. Mean Magnitude 

In Fig. [4] we plot the mean V magnitude of the system 
as a function of time, as measured by AS AS. The empty 
squares represent the weighted mean of the observations, 
taken over four months in each observing season and the 



error bars represent the square r oot of the variance of the 
weighted mean, computed as (see lFrodesen et aD . ll979l Eq. 
9.12): 



Var — 
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Conjunction epocli [HJD] 
Period [days] 
asini [R0] 

Systemic velocity [kms^^] 
Primary RV amplitude [kms~^] 
Secondary RV amplitude [kms~^] 
Orbital excentricity 
Misin^ j [M0] 
Mi[M0] 
M2sin^j [Mq] 
M2 [Mq] 



2452779.422 ± 0.005 
1.87844246 ± 0.0000049 
5.665 ± 0.035 
18.81 ±0.32 
74.58 ±0.68 
78.07 ±0.66 
0.0072 ± 0.0070 
0.3541 ± 0.0068 

> 0.3646^ 
0.3382 ± 0.0068 

> 0.3482 " 



" Infered from the absence of eclipses in the lightcurve 
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Fig. 5. Lomb-Scargle periodogram of the ASAS data. The 
main peak occurs at Pphot — 1.876644 days, and has a False 
Alarm Probability of about 0.1%. The horizontal dashed 
lines show the values of the False Alarm Probability. 
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Fig. 4. Mean V magnitude of GL 375 as measured by the 
All Sky Automated Survey as a function of time. The mean 
was taken every four months in each observing season. 



where ai is the photometric error of ith point, which is pro- 
vided by ASAS, and n is the total number of observations 
analizcd. Note the lack of observations during part of the 
2002 season, as well as in the entire 2005 and the beginning 
of 2006 seasons. 

A maximum is discernable around the year 2002.5 
{HJD w 2452430), and a minimum in 2004 {HJD « 
2453050). A few additional minima seem to be present near 
2002 {HJD « 2452200) and 2003.25 {HJD « 2452700). 
Unfortunately, the lack of observations during the last part 
of 2004 and 2005 precludes any further analysis at this 
point, since a maximum could have occurred in 2005. We 
return to this issue in Sect. 14.31 

4.2. Rotational Modulation 

We calculated the Lomb-Scar gle periodogram (jScargld 
119821: iHorne fc Bahunad . I1986D for the ASAS data. This 
periodogram is a method to estimate the power spectrum 
(in the frequency domain) , when observation times are un- 
evenly spaced, and normalizes the spectrum in such a way 
that it is possible to estimate the significance of the peaks. 
The amplitude of the periodogram at each frequency point 
is identical to the equation that would be obtained esti- 
mating the harmonic content of a data set, at a given fre- 
quency, by li near least-squares fitting to an harmonic func- 
tion of time ([Press et al.l . ll993 ). We show this periodogram 
in FigO 



As can be seen, there is a distinct peak corresponding to 
a period of Pphot = 1.876644 days. For this peak, the False 
Alarm Probability is about 0.1%. The excellent False Alarm 
Probability indicates the presence of a harmonic variability 
within the data set with a period resembling very closely 
the measured orbital period. 

We believe this harmonic variability is produced by 
spots and active regions in the stellar surface carried along 
with rotation. This would imply that the rotational and 
orbital periods are synchronised, as is expected for such a 
close binary. 

To verify if this is indeed the case, we phased the data 
to the obtained period, for two different seasons (Fig. [6|). 
The sinusoidal shape of the variation is evident, although 
the amplitude of the modulation is different in both cases: 
between 2002.5 and 2003.5 (upper panel) the amplitude is 
about 30 mmag, and falls to about 15 mmag during 2006. 
This could indicate a different area covered by starspots or 
active regions, and therefore different activity levels, in each 
epoch. Note that during the period of minimum modulation 
the mean brightness is smaller by about 20 mmag. 

Note also that whether the period of maximum rota- 
tional modulation corresponds to an epoch of maximum 
or minimum activity depends on the filling factor. On the 
other hand, the variations in the mean brightness can be 
dominated by the active regions (as in the Sun) or by the 
spots (as in younger stars), which would cause the stars 
to be either brighter or fainter during periods of maxi- 
mum activity, respectively. Therefore, considering that dur- 
ing the period of minimum modulation the obseved system 
is fainter (see Fig. [6]), we are left with two possibilities. 
Either the system is dominated by active regions and the 
filling factor during maximum activity is below 50% (as in 
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Fig. 6. ASAS photometry phased to the orbital period for two different epochs. Left: 2002.5 - 2003.5 Right: 2006 
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Fig. 7. Variation in the binned photometric data as mea- 
sured from a sinusoidal fit (see text) as a function of time. 
A peak is clearly visible near the year 2003. 

the Sun) , or is the spots that control the photometric vari- 
ations and the maximum filling factor is over 50%. 

We found no evidence of eclipses in the photometry, 
and since the ASAS data have complete phase coverage in 
many periods, this could not be an eclipsing system, and 
we can therefore set an upper limit to the inclination angle. 
Assuming a stellar radius of 0.4 R0 for both components, 
the inclination angle should be larger than 82°. In turn, this 
implies lower limits for the stellar masses. Mi > 0.3646 Mq 
and A/2 > 0.3482 Mq. Note that the assumed spectral type 
of the stars remain consistent with these values. 

To study the evolution of the activity level, we studied 
in greater detail the changes in the amplitude of the pho- 
tometric modulation as a function of time. Assuming that 
the modulation is produced by active regions and spots in 
the stellar surface, its amplitude should trace the amount 
of inhomogeneities in the surface, and therefore the level of 
magnetic activity. 

We binned the ASAS measurements every four months 
and phased the data on each bin to the orbital period. Bins 
in which less than 7 observations were present were dis- 
carded. The phased data was then fitted to a sinusoidal 
model of the form: 



Oo + CLs sin(27riy9) -|- cos(27r(p) 



(1) 



where ip is the orbital phase. The resulting fit parameters 
and their covariance matrix were employed to measure the 
amplitude A of the variation and its standard error a a in 
every time bin: 



A ~ \/ fls^ + 
a A = A^^[veLT{as) + ySiT{ac) + 2 as ttc 



1I/2 



(2) 
(3) 



where var(aa;) denotes the variance of the parameter 
and cov is the covariance of the parameters. The results 
are shown in Figure [7l where the time assigned to each 
measurement was the mean time of the photometric data 
considered. 

It is interesting to compare the present Figure with 
Fig. [3) The rotational modulation also exhibits a maximum 
around the year 2002.5, as in the previous section, and there 
are also two marked minima around the year 2002 and 2004, 
although the minima of amplitude of the rotational mod- 
ulation are more significant. Consider that although these 
measurements originate in the same dataset, the presence 
of minima and maxima at the same time is meaningful, 
since the measurements are of different nature. 
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Fig. 8. Comparison of the fluxes in the Call K line for both 
components. The error bars correspond to a 10% error in 
the line fluxes, and the dotted line is the identity relation. 



4.3. Magnetic activity 

To study the magnetic activity of the system we use the 
fiux in the Call K line, a well-known activity indicator. We 
did not use the Call H line, since it is contaminated by the 
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Table 3. Flux in the Ca II K line for both components of 
the system 
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Fig. 9. Line fluxes in Ha and in the Ca II K line for the 
primary (empty squares) and the secondary (filled circles) 
components of the system. It can be seen that the line fluxes 
are uncorrelated for both stars. The error bars correspond 
to a 10% error in the line fluxes. 



nearby Balmer line He, which we expect to be in emission 
in these active stars. 

Since we have separated the spectra from both compo- 
nents, we should be able to distinguish whether the changes 
observed in the photometric measurements belong to one 
of the components or to both. The spectroscopic measure- 
ments also provide further information about the active 
regions present in the stellar surfaces. For example, com- 
bining the photometric data with measurements of activity 
proxies it should be possible to discern whether active re- 
gions or stellar spots dominate the photometric variations, 
and whether the filling factor is larger or smaller than 50%. 
Furthermore, since we have observations for the year 2005, 
a further detailed study of the system variability is possible. 

We carefully examined the separate spectra and dis- 
carded those where reduction problems were present or 
where the S/N was too poor. We also discarded spectra 
which exhibited signs of a transient event, such as a flare. 
To do this we employed the individual observations (see 
Sect. [2]) and compared the amplitudes of the K line fluxes in 
one individual spectrum with those in the other. If no flare 
were present, no significant difference should be observed, 
since the time between observations is about an hour. On 
the other hand, if a strong transient event is present during 
our observations, we should observe differences in the K line 
fiuxes from one spectrum to the next. In the last column of 
Table[T]the discarded observations are marked with a ! sign. 
In some cases, generally due to changing weather conditions 
or the presence of flares, only one of the individual spectra 
had to be discarded. In those letter next to the sign 

indicates which spectra was discarded, either a or b, for the 
first and second one, respectively. Additionally, the spec- 
tra taken during the night of March 2nd 1999 (marked §) 



were discarded because the spectral range does not include 
the Ca II K line, due to the use of a different instrumental 
setup during that night. The flux measurements were car- 
ried out in the remaining individual spectrum, after check- 
ing that the regions of interest were free of cosmic rays. 
Unfortunately, one of the discarded spectra belong to the 
year 2005 and another one to the end of 2004. 

The fluxes in the K line were computed in a 1 A win- 
dow using the IRAF task SBANDS, which measures the 
flux by summing each pixel in the deflned interval. At the 
edges of the bandpass it considers only the fraction of the 
pixel laying within the bandpass. The results are presented 
in Table [31 In Figure [5] we plot the K line flux for both 
components of the system. The error in the line fluxes was 
conservatively estimated as 10%. It can be seen that the 
fluxes are strongly correlated (r = 0.95) and, therefore, the 
activity variations of the components should be in phase. 
For this to occur, the magnetospheres of the stars must be 
strongly interacting. Vahia (1995) simulated the interaction 
of the magnetic fleld of binary stars. In that work, the stars 
are modelled as perfect dipoles surrounded by vacuum and 
it is shown that it should be quite common to find one 
of the components located completely inside the magne- 
tosphere of the other. Various works describe mechanisms 
by which magnetic activity can be enh anced in interacting 
systems like this one (see, for example, IZaaarashvili et all 
2002) . This kind of mechanisms could in principle explain 
the high level of emission present in the lines used as activ- 
ity proxies, such as the Ca II K line. In Figure [5] the dotted 
line is the identity relation. Note that the line fiuxes of the 
primary component are systematically larger than those of 
the secondary. The mean flux of the primary is around 13% 
larger than on the secondary, and the difference appear to 
be greater for larger fluxes. 

In Fig. [To] we plot the flux in the K line as a function of 
time. The flUed (empty) diamonds and solid (dashed) lines 
represent observations of the primary (secondary) compo- 
nent. Distinct maxima are present for both components 
around the years 2002 and 2004. There are also significant 
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2000 2002 ^^""^ 2004 2006 




Fig. 10. Flux in the Call K line as a function of time for 
the primary (solid line) and secondary (dashed line) com- 
ponents of the system. The error bars correspond to an 10% 
error in the line fluxes. 



minima at 2001.5 and 2003.5, and a less clear one at around 
2006. Apparently, the periods with little activity (minima) 
are much shorter than those with higher level of chromo- 
spheric emission, and are separated by about 2 years. 

Another activity indicator usually employed in the 
study of M stars is the flux in Ha. The flux in this line 
is generally considered to correlate well with the flux in 
the Call H and K lines and has the advantage of being lo- 
cated in a redder wavelength range than the Ca II lines, 
where M dw arfs are brighter. However , it has been recently 
reported bv lCincunegui et al.l (|2007br ) that the correlation 
between Ca II and Ha is not always valid. In that paper we 
found that while some stars exhibit correlations between 
Ha and the Ca II lines, the slopes change from star to star. 
Furthermore, other stars exhibit anticorrelations and there 
are also cases where no correlation is present. Therefore, we 
checked whether a correlation existed for the components 
of Gl 375. In Fig. [9] we plot the lines fluxes in Ha vs. those 
in the Ca II K line for both components of Gl 375. It can 
be seen that within the errors, none of the components ex- 
hibit any sign of correlation between the lines fluxes and, 
therefore, this line cannot be used as an activity indicator 
for these stars. 

In the next Section we study in detail if the system 
presents a definite activity cycle, like the one present in the 
Sun. To do this, we compare the three observed magnitudes. 

4.4. Periodicity 

In Fig. [TT]we plot together the three measured quantities: 
mean magnitude (empty squares and solid line) , amplitude 
of the rotational modulation (filled square and dotted line), 
and the K line flux (empty circles and dashed line). Since 
the fluxes in the K lines of both components are in phase, 
for clarity in the graph we plot the total flux of the system. 
Also to avoid crowding the graph, we have not included the 
error bars. The K line flux curve has been shifted by 140 
days to coincide with the other curves. The arrow in Fig. 1111 
indicates this displacement. 

This timelag between photometric and magnetic vari- 
ations has already been observed for stars of differ- 
ent spectral types, from /3Com (GOV) to eEri (K2V), 
including binary syste ms like ^Boo (see , for example, 
iGray &: I3aliunas.. ,1995; ,Grav et al.l . Il996bl lal). In particu- 
lar, Grav et al.l ( 199631) showed that the time lag between 
magnetic variations, and variations observed in photometry 



Fig. 12. Total flux of the system in the Call K line as a 
function of time. The error bars correspond to a 10% error 
in the line fluxes, and the solid and dashed lines are the 
best sinusoidal fit and the ±3cr deviations, respectively. The 
period obtained from the fit is 763 days. 



and temperature is anticorrelated with effective tempera- 
ture (see their figure 8). For the coolest star in their sample 
(eEri) the time lag of temperature variations is about 0.3 
years. The timelag present in the Gl 375 system is compa- 
rable to this value. 

As mentioned before, the position of the maxima and 
minima of the mean magnitude (empty squares and solid 
line) and those observed in the amplitude of the rotational 
modulation (filled square and dotted line) coincide very 
well. It can also be seen that, once displaced, the behaviour 
of the K line flux also agrees very well with the other two ob- 
servables. The maximum around the year 2002.5 is present 
for all curves, as well as the minima near 2001.8 and 2004. 
Note also that the distinct minimum in the line flux present 
around the year 2006.2 is also present in the curves of the 
other two observables, although the minima are much less 
significant, and that the slight decrease in magnitude and 
amplitude around 2003.2 coincides with a decrease in the 
line flux around 2002.8. Unfortunately, no spectrum was 
obtained between the one at 2002.8 and the one at 2003.5. 

The agreement between the behaviour of the three ob- 
servables is remarkable, because of the different nature of 
the observations, and the different instruments and sites in 
which they were obtained. Although no clear periodicity 
can be concluded from the present data set, the main min- 
ima present in all three curves appear to occur periodically 
every 2.2 years. This hints to the presence of a periodic 
activity cycle in this system, although further observations 
are necessary. 

To further study this point we fit the period, ampli- 
tude and phase of a sinusoidal function to the total flux us- 
ing a modification of the Levenberg-Marquardt algorithm 
available in the SciPy library for Python. The fit is shown 
in Fig. [T^l together with the ±3 a deviations, for the con- 
servative estimate of 10% in the line fluxes. The obtained 
period is P = 763 days, with a formal error of less than 
5 days. However, the fit is extremely sensitive to some of 
the points. For example, if the lowest point, located near 
HJD - 2450000 = 3700 is not considered in the fit, the ob- 
tained period changes to 786 days. Therefore, we consider 
that a more realistic error for the period would be of 5%, 
i.e. about 40 days. 
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Fig. 11. Mean magnitude (black solid line and empty squares), amplitude of the rotational modulation (red dotted line 
and filled squares) and flux in the Call K line (blue dashed line and empty circles) as a function of time. The arrow 
indicates the temporal displacement of 140 days that have been applied to the line flux (see text). 



5. Conclusions 

We studied the spectroscopic binary system Gl 375 us- 
ing medium resolution echelle spectra obtained at the 
Argentinian observatory CASLEO since the year 2000, and 
photometric measurements in the V band from the ASAS 
database. We found the system exhibits enhanced levels of 
chromospheric emission, as is expected for binary systems. 

From each observed spectrum we obtained the individ- 
ual spectra and RV of each component. Wc found that the 
two components are dM3.5 stars, and that they orbit each 
other every ~ 1.88 days in a circular orbit, with a minimum 
radius asinz 5.7 Kq. The minimum masses obtained (see 
Table [2]) are similar to the expected values for this type of 
stars and, therefore, the orbital inclination angle is expected 
to be ~ 90°. However, the system does not exhibit eclipses 
and, therefore, we are able to set an upper limit to the incli- 
nation angle (i < 82°), which in turn implies a lower limit 
for the masses of the components (Mi > 0.3646 Mq and 
M2 > 0.3482 Mo). 

As expected for stars in such a close orbit, their rota- 
tional periods are synchronised to the orbital period. This 
is observed as a modulation in the V-band light curve with 
the same period as the orbit. Additionally, the amplitude 
of this modulation is clearly seen to change with time in 
longer timescales, which we believe is due to differences in 
the area of the stars covered by the active regions. 



Therefore, by studying the evolution of the amplitude 
of the rotational modulation we found that the system ex- 
hibits a roughly periodic behaviour of 2.2 years (or ^ 800 
days). The same period was found in the mean magnitude of 
the system and in the flux of the Ca II K line, although the 
Ca flux variations occur 140 days ahead of the photomet- 
ric ones. This behaviour was observed previously in other 
single stars and binary systems, for which the lag between 
temperature variations and line fluxes was reported to an - 
ticorrelate with effective temperature (|Grav et al.l . Il996al) . 
The measured lag for this system fits roughly in that trend, 
indicating that the same process may be responsible for the 
phase shift between observables. 

The agreement between the behaviour of the three ob- 
servables is remarkable, because of the different nature of 
the observations, and the different instruments and sites 
in which they were obtained. Additionally, we performed 
a non linear fit to the total Ca II K fiux using a sinu- 
soidal function, and found the least-square estimation for 
the period is P = 763 days. If confirmed, this would be 
the second detection of an activity cycle in an M star, the 
first being a recently reported cycle o f aroun d 450 days 
in Proxima Centauri (jCincunegui et all . l2007a[ ) . It is inter- 
esting to compare briefly both cases in relation with the 
2.5-years lower li mit for the pe riod of activity cycles in so- 
lar stars given bv lBaliunas et a l. (1995). Although slightly 
shorter, the period of the Gl 375 system is consistent with 
this limit, which is reasonable, since Gl 375 A and Gl 375 
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B are dM3.5 stars, and should therefore have a radiative 
core surrounded by a convective outer envelope. Therefore, 
in the Gl 375 system the standard dynamo could in princi- 
ple explain the cyclic behaviou r observed, as in the e arlier 
stars included in the study bv iBaliunas et alj (|1995[ ). On 
the other hand, the period of Proxima Centauri is well be- 
low the 2.5-years cutoff. Since Proxima Centauri is a dM5.5 
star, and should therefore be completely convective, a pos- 
sible interpretation of this fact is that different mechanisms 
are responsible for the production of the activity cycle in 
Gl 375 and in Proxima Centauri. Further observations of 
these systems, as well as the discovery of activity cycles in 
other M stars should provide a solid observational base to 
develop a model capable of explaining the observations. 

Another interesting result of this work is that the ac- 
tivity of Gl 375 A and Gl 375 B, as measured in the flux 
of the Ca II K lines, are in phase. The levels of chromo- 
spheric emission of both components exhibit an excellent 
correlation, which implies that a magnetic connection exists 
between both components, as is shown in numerical simu- 
lations by other authors. Due to its vicinity and relative 
brightness this system presents an interesting opportunity 
to further study this type of interaction. 

Finally, since the chromospheric emission is larger when 
the system is brighter, and this coincides with period of 
maximum amplitude of the modulation, we conclude that 
the photometry of these stars is dominated by bright active 
regions rather than dark spots and that the filling factor is 
below 50%, very much like in the Sun. 
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